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I review our theoretical understanding of thermonuclear flashes on accreting neutron stars, concentrating on 
comparisons to observations. Sequences of regular Type I X-ray bursts from GS 1826-24 and 4U 1820-30 are very 
well described by the theory. I discuss recent work which attempts to use the observed burst properties in these 
sources to constrain the composition of the accreted material. For GS 1826-24, variations in ot with accretion rate 
indicate that the accreted material has solar metallicity; for 4U 1820-30, future observations should constrain the 
hydrogen fraction, testing evolutionary models. I briefly discuss the global bursting behavior of burst sources, 
which continues to be a major puzzle. Finally, I turn to superbursts, which naturally fit into the picture as 
unstable carbon ignition in a thick layer of heavy elements. I present new time- dependent models of the cooling 
tails of superbursts, and discuss the various interactions between superbursts and normal Type I bursts, and what 
can be learned from them. 



1. INTRODUCTION 

Type I X-ray bursts were discovered in the 
1970's, and quickly understood as being due to 
unstable burning of hydrogen and helium on the 
surface of an accreting neutron star (for reviews 
see |34I48| '1. Fuel accumulates for hours to days, 
and then ignites and burns in a « 10-100 s burst 
with typical energy w 10'^^-10''° ergs. 

In recent years, there has been a lot of ex- 
citement in the field. First, the Rossi X- 
Ray Timing Explorer (RXTE) discovered nearly- 
coherent oscillations during bursts, with frequen- 
cies in the range 270-620 Hz gS]- Believed 
to be rotationally-modulated brightness asymme- 
tries on the neutron star surface, the oscillation 
frequency tells us the neutron star spin, and po- 
tentially about the spreading of burning, and the 
surface magnetic field 0. 

Second, long term monitoring of bursters with 
BeppoSAX/WFC (Verbunt, this volume) and 
RXTE/ ASM and PC A has revealed a new class 
of rare, extremely energetic, and long duration 
bursts now known as "superbursts" ([301; Kuulk- 
ers, this volume), as well as a new class of low 
luminosity burst sources discovered through their 
burst rather than persistent emission Cor- 
nelisse, this volume). Particularly important has 



been the large amount of data accumulated on 
bursters by the Wide Field Camera (WFC) sur- 
vey of the Galactic center, which has also enlarged 
our sample of normal Type I bursts 14 . 

In this article, I review our theoretical under- 
standing of thermonuclear flashes on accreting 
neutron stars, concentrating on how well it com- 
pares to observations, and what we can learn from 
such a comparison. 

2. HOW THE ACCRETED HYDROGEN 
AND HELIUM BURNS 

The physics of nuclear burning on the surface 
of an accreting neutron star has been discussed 
by many authors, comprehensive reviews can be 
found in refs. |34I4I48| . Here I briefly review how 
the accreted hydrogen and helium burns, and the 
expected burst properties as a function of accre- 
tion rate. In particular, I concentrate on the dif- 
ferences between hydrogen and helium burning, 
which are crucial for understanding the regimes 
of burning, and for comparing to observations. 

The timescale for burning — whereas he- 
lium burning proceeds rapidly via strong interac- 
tions, hydrogen burning is relatively slow, since 
beta decays are required to convert protons to 
neutrons. This leads to a saturation of the CNO 
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burning rate at high temperatures T >8 x 10'' K 
(the so-called "hot CNO cycle"; [23 ), when pro- 
ton captures in the cycle occur more rapidly than 
the subsequent beta-decays. At high accretion 
rates, hot CNO burning during accumulation of 
fuel depletes hydrogen prior to the runaway. The 
time to burn all the hydrogen in a given fluid el- 
ement is 

'-"-(^)(ff). 

where Z is the CNO mass fraction. 

The amount of hydrogen vs. helium present 
at ignition of the burst has a direct impact on 
the burst lightcurve. If the helium fraction is 
large, the rapid energy release at the beginning of 
the burst creates a luminosity exceeding the Ed- 
dington luminosity, driving outwards expansion 
of the photosphere ( "photospheric radius expan- 
sion") 123. If hydrogen dominates the compo- 
sition however, it burns after ignition via the rp- 
process |57I42| . Starting with seed nuclei made by 
helium burning, this comprises a series of proton 
captures and beta-decays involving nuclei close to 
the proton drip line and extending well beyond 
the iron group. This can significantly delay the 
cooling of the burning layer, giving an extended 
burst tail, and burst durations ~ 100 s |27l60j . 

Thermal stability — The nature of the insta- 
bility is a thin shell flash |2H1, in which tempera- 
ture perturbations drive runaway heating through 
temperature sensitive nuclear reactions. Hot 
CNO hydrogen burning is thermally stable, since 
the burning rate does not respond to tempera- 
ture fluctuations. However, at the lowest accre- 
tion rates M < 0.01 MEdd\ T < 8 x lO'^ K in 
the accumulating layer and the usual CNO cy- 
cle operates, allowing unstable hydrogen ignition 
and leading to "hydrogen-triggered" bursts. For 
M > 0.01 A/Edd, the hydrogen burns stably, and 
helium burning at a density > 10^ g cm^'^ drives 
the thermal instability. 

Helium triggered bursts are of two types, de- 
pending on the time to accumulate the critical 

^depending on metallicity, see |1] for the scalings. I take 
A^Edd = 1-7 X 10~* Mq yr~l, the Eddington accretion 
rate for a 10 km neutron star accreting solar composition 
material. 



column of fuel (or equivalently the recurrence 
time of the bursts irocur)- If ^rccur > iff J hydrogen 
burns away and ignition occurs in a pure helium 
layer below a steady burning hydrogen shell, giv- 
ing a "pure He" burst; if trecur < tn, hydrogen 
is present throughout the layer when the helium 
ignites, giving a "mixed H/He" burst. 

For M > MEdd , helium burns at a temperature 
T > 5 X 10® K for which the temperature sensi- 
tivity of the triple alpha reaction is less than the 
temperature sensitivity of the cooling. The burn- 
ing is then thermally stable, and the fuel burns at 
the rate it accretes. The accretion rate at which 
burning stabilizes and its dependence on metal- 
licity is still not well-determined theoretically. 

Nuclear energy release — hydrogen burn- 
ing releases substantially more energy per gram, 
Qnuc ~ 7 MeV per nucleon compared to only 
w 1.6 MeV per nucleon from helium burning. Ob- 
servationally, the energy release per gram is mea- 
sured by the parameter a, the ratio of persistent 
fluence between bursts to burst fluence, 

^ CF.dt ^ (GM/R) 

CF^dt^ Qnuc ' 

where Fp is the persistent flux from accretion, .Fb 
is the burst flux. At is the interval from the begin- 
ning of one burst to the next, and GM/R « 200 
MeV per nucleon is the gravitational energy re- 
lease. We expect a « 40 if the material burned 
in the burst is solar composition, compared to 
a >; 100 for pure helium. (These numbers are 
somewhat uncertain given that not all of the 
accretion energy necessarily appears in X-rays, 
and the X-ray emission is potentially anisotropic). 
Measurement of a values in the expected range 
was an important initial confirmation of the ther- 
monuclear flash model for Type I bursts. 

The large energy release from hydrogen burn- 
ing plays an important role in setting the condi- 
tions for ignition of the flash. Figure^shows illus- 
trative temperature profiles as a function of col- 
umn depth (g cm~^) into the star at the moment 
of ignition for each of the four different burn- 
ing regimes for the case of solar metallicity. The 
dashed and dotted curves are ignition curves for 
hydrogen and helium, calculated using a local ap- 
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proximation for the cooling rate following ref. JHI 
(I refer to that paper for details, also see [37] for 
recent calculations of the global eigenfunctions) . 

In the mixed H/He burst regime for so- 
lar metallicity, the ignition conditions are well- 
determined, since the temperature profile is set 
by hot CNO burning within the accumulating 
layer. The ignition column depth is « (1- 
2) X 10'* g cm^^, roughly independent of accretion 
rate, so that bursts in this regime are expected to 
have w constant energy, and a recurrence time 
which scales close to cx 1/M. 

For low metallicity mixed H/He bursts, pure He 
bursts, and hydrogen-triggered bursts, the inter- 
nal heating from hot CNO burning is either much 
less or absent altogether. The temperature pro- 
file is then set by compressional heating or heat- 
ing from below, either heat from reactions occur- 
ring deep in the neutron star crust 6 , or residual 
heat from or burning of leftover fuel in the ashes 
of a previous burst (51I60| . The ignition con- 
ditions are therefore sensitive to accretion rate, 
and might be expected to show variations from 
burst to burst (i.e. less regular bursting). This 
is particularly true for pure He bursts because of 
the very shallow slope of the ignition temperature 
with column depth. 

3. COMPARISON WITH OBSERVA- 
TIONS 

3.1. Regular bursting I: GS 1826-24, The 
"Clocked Burster" 

The transient source GS 1826-24 was discov- 
ered by Ginga, and regularly monitored by Bep- 
poSAX, which revealed it to be an extremely reg- 
ular Type I burster. Ubertini et al. jSSj, who 
dubbed this source the "clocked" burster, found 
the burst recurrence time was close to 6 hours 
with a dipersion of only « 6 minutes. Further 
monitoring showed that the recurrence time de- 
creased as the source brightened P^ . 

Based on the energetics and recurrence times, 
as well as the long duration of the bursts, Bild- 
sten proposed that this was a textbook case 
of mixed H/He bursting. This is the expected 
burning regime at the inferred accretion rate of 
0.1 AfEdd and observed recurrence times of < 6 
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Figure 1. Illustrative temperature profiles (solid 
curves) in the different burning regimes, for solar 
metallicity. The dotted and short-dashed lines 
are helium ignition curves (helium mass fraction 
Y = 0.3 and F = 1.0 respectively); the long- 
dashed line is the hydrogen ignition curve. 



hours, and is consistent with the measured a = 
40-60. In addition, the inferred ignition mass is 
as expected from theory. The burst energy of 
« 5 X 10'^^ erg implies an ignition mass of w 10^^ g 
for a nuclear energy release 5 x 10^^ erg g~^ 
(equivalently, this is the mass accreted in a recur- 
rence time at the observed M). The correspond- 
ing column depth is j/ign ~ lO'* g cm^^, exactly 
as expected for mixed H/He ignition with solar 
metallicity {Z = 0.02) (Figure □ see also Table 2 

of [El). 

Given the excellent agreement with theory, a 
natural question is whether we can use the bursts 
to constrain properties of the source [23. For 
mixed H/He bursts, the ignition depth is mainly 
sensitive to the metallicity, because the temper- 
ature profile of the accumulating layer is set by 
hot CNO burning. This suggests an argument to 
constrain the metallicity: if it was less than so- 
lar, the reduced CNO heating would give a larger 
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Figure 2. Ratio of persistent to burst luminosity 
a = Lp/Lb (equation |2Jl, calculated from RXTE 
observations between 1997 and 2002. Error bars 
represent the estimated la uncertainties. From 
Galloway et al. (2003) (ref.[211). The changing 
a with M indicates that fuel is burning between 
bursts, as expected for hot CNO hydrogen burn- 
ing with solar metallicity. 



ignition column depth, and a larger recurrence 
time and burst energy than observed. However, 
the uncertainty in the relation between M and 
the X-ray luminosity Lx allows lower metallicity 
models to be accommodated by adopting a larger 
M (the larger column is then accreted in the ob- 
served recurrence timc)^. Additionally, the time- 
dependent simulations of Woosley et aLlS]] show 
that at low metallicity, a new heat source takes 
over. Production of CNO elements in the tail 
of the previous X-ray burst leads to leftover hy- 
drogen burning beneath the newly accreted layer 
("compositional inertia"), reducing the sensitiv- 
ity to metallicity (see also 127!). 

Determining the metallicity requires an addi- 
tional constraint on the models. Galloway et 
al. recently concluded that the metallicity in 

^In fact, the best approach is to turn this around: given 
a metaUicity, we "measure" the local accretion rate onto 
the star by matching the observed recurrence time. 
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Figure 3. The average observed lightcurve from 
GS 1826-24 (dotted line; ref. EH) compared with 
theoretical simulations (solid and dashed lines; 
models ZM and zM of ref. (HOI)- The observed and 
theoretical recurrence times are both « 4 hours. 



the accreted layer was solar by looking at vari- 
ations of burst properties with M . They present 
data for 24 bursts observed by RXTE between 
1997 and 2002. The a value was found to decrease 
by ~ 10% as accretion rate increased by ~ 50% 
(Figure This is exactly as expected from 

models with solar metallicity [TB], since as ac- 
cretion rate increases and recurrence time drops, 
less time is available to burn hydrogen between 
bursts. Models with metallicity much less than 
solar are hard to reconcile with the a variations. 

Another constraint comes from the burst 
lightcurves, which show a long w 10 s rise and 
a long « 100 s tail^. Theoretical models of 
mixed H/He bursts (e.g. |27l43j ) show long burst 
tails, leading Bildsten to argue that an ac- 
tive rp-process powers the bursts from GS 1826- 
24. Figure|2|compares the average observed burst 
lightcurve with the recent simulations of Woosley 
et al. These simulations are the first to in- 

corporate a large nuclear reaction network that 
can follow the detailed rp process nucleosynthe- 
sis together with resolved vertical structure. The 

^The burst lightcurves were remarkably similar during this 
period, with only a small variation as M increased with 
time. This, together with the extremely regular nature 
of the burning, argues for complete burning of the fuel in 
each burst, and complete covering of the stellar surface. 
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agreement between the observed lightcurve and 
solar metallicity model is remarkable, support- 
ing the findings of Galloway et al. 24 . A more 
detailed comparison with the simulations is in 
progress. The predicted burst lightcurves are sen- 
sitive to the nuclear physics input |60| . and so 
there is potentially much to learn. 

The ignition column for mixed H/He bursts 
is predicted to be roughly independent of accre- 
tion rate, since hot CNO burning sets the tem- 
perature in the accreted layer (§2). Galloway et 
al. [22 find that i^ccur oc A/-i 05±o o2^ implying 
a slight decrease in the ignition column with M, 
yign oc M~" '^^. This remains to be understood; a 
simple argument suggests a slight increase should 
be seen, since less helium is produced by hot CNO 
burning for a shorter recurrence time, requiring a 
slightly thicker layer to ignite. Two possible ex- 
planations for the discrepancy are (i) extra heat- 
ing, e.g. from thermal or compositional inertia ef- 
fects |51I60| , or (ii) the fraction of the star covered 
by fuel changes with M . Future comparisons with 
time-dependent simulations or improved spectral 
models should test these ideas. 

3.2. Regular bursting II: 4U 1820-30 

The ultracompact binary 4U 1820-30 (orbital 
period 11.4 mins '47') is an interesting source in 
which to study Type I bursting behaviour. It has 
a known distance, being located in the metal-rich 
globular cluster NGC 6624 ([Fe/H] « -0.4, dis- 
tance 7. 6±0. 4 kpc) It undergoes a character- 
istic 176 day accretion cycle ^U] during which the 
accretion rate varies by a factor of 3 or more. Reg- 
ular sequences of Type I bursts are observed dur- 
ing the low state, with recurrence times of ~ 2- 
4 hours |11I12I26I14) . In the high state, bursts 
become irregular or disappear altogether 10 14 . 
This source also showed the most energetic and 
luminous superburst I49j . 

The extremely short orbital period implies a 
hydrogen-deficient companion, and observations 
of bursts support this picture, with a « 120 
|26| as expected for helium rich material. Dif- 
ferent evolutionary models predict that the ac- 
creted matter is either pure helium j55lll41j . or 
contains a small amount of hydrogen (perhaps 5- 
35% by mass; [SaHSI)- Recently, I looked at 
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Figure 4. Predicted burst energy and recurrence 
time as a function of assumed local accretion 
rate m for 4U 1820-30 (ref.jljj). Choosing m al- 
lows both pure helium models (dotted and short- 
dashed lines) and models with hydrogen (long- 
dashed line 10% hydrogen; solid line 20% hydro- 
gen) to fit the observed values (indicated by error 
bars and placed at the value of m inferred from 
the X-ray luminosity). 



the question of whether the composition of the 
accreted material (specifically the hydrogen frac- 
tion) might be inferred from the burst proper- 
ties |15|. This is important to test intermediate- 
mass binary evolution models for X-ray binaries 
(e.g. |39p which predict a small amount of hydro- 
gen may be present. 

The inferred accretion rate when bursts are 
seen from 4U 1820-30 is M w 2.4 x 10"^ Mq yr'^, 
giving an ignition column depth « 2 x 10^ g cm~^. 
Bildsten 3 found that this compares well with 
time-dependent models of pure helium burning, 
although for a slightly hotter base temperature 
than expected. The reason for this can be seen by 
inspecting Figure Q Ignition at y w lO'* g cm~^ 
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requires heating equivalent to hot CNO burning 
with Z ^ 0.01. In the pure hehum case, this heat- 
ing must be provided by a heat flux into the layer 
from below, and translates to a flux equivalent to 
K, 0.4 MeV per accreted nucleon, higher than the 
expected crust luminosity 

If hydrogen is present with mass fraction X 
0.1, hot CNO burning during accumulation could 
provide the required heating ^21 (this requires 
Z > 3 X 10-3, hkely satisfied for 4U 1820-30). 
This is illustrated in Figure 0] which shows the 
predicted recurrence times and burst energies for 
models with and without hydrogen, compared to 
observed values. Models with hydrogen do much 
better for the accretion rate inferred from the X- 
ray luminosity. However, it is not possible to 
reach a firm conclusion because increasing the ac- 
cretion rate by only a factor of two allows pure 
helium models to come into agreement. 

As for GS 1826-24, an additional constraint 
is required to determine the accreted composi- 
tion. One possibility is to look for variations of 
the burst energetics with accretion rate. For 10% 
accreted hydrogen, the time to burn all the hy- 
drogen is « 3 hours, comparable to the recur- 
rence time. This leads to sa 10% variations in 
burst fluence as recurrence time varies. This is 
a promising way to constrain the accreted com- 
position, but will require further observations of 
sequences of regular bursting with different recur- 
rence times. Another possible constraint comes 
from simultaneous modelling of Type I bursts 
and superbursts, I discuss this further in H4.3I 
It would also be interesting to compare burst 
lightcurves with time-dependent simulations. 

3.3. Global bursting behaviour 

The properties of regular burst sequences from 
both GS 1826-24 and 4U 1820-30 agree very well 
with theoretical expectations, and we have dis- 
cussed to what extent we can learn about the 
accreted composition in these cases. However, 
it has been pointed out many times that the 
global bursting behavior is not understood, and 
for many bursters is opposite to that predicted 
by theory |56I5| . EXOSAT observations of sev- 
eral Atoll sources showed that as X-ray luminos- 
ity increased, burst properties changed from reg- 



ular, frequent bursts (trccur ~ hours) with long 
durations (« 30 s) and low a w 40, to irregu- 
lar, infrequent (irccm- > 1 day), short (< 10 s) 
bursts with high a > 100 and often showing 
photospheric radius expansion (PRE) (e.g. 56.)- 
RXTE and BeppoSAX observations have con- 
firmed this result, with particularly good cover- 
age for the transient source KS 1731-260 |35ll4j . 
Cornelisse et al. found that observations of 
nine bursters were consistent with this pattern of 
bursting, with a universal transition luminosity 
Lx w 2 X 10^^ erg s"^ (although not all showed 
the transition from one burst type to another). 

The short bursts are naturally interpreted as 
pure He bursts; the long bursts as mixed H/He 
bursts*. However, theory predicts a transition 
from pure He bursts to mixed H/He bursts as 
accretion rate increases rather than the other 
way around (§2). Several explanations have been 
proposed, all of which require further investiga- 
tion. In the context of ID models, mixing of 
fuel by Raylcigh- Taylor [S^ or shear instabilities 
|22j might allow more rapid burning of hydrogen; 
Woosley et al. have again emphasised the im- 
portance of "compositional inertia" , that bursts 
ignite on the ashes of previous bursts; Narayan & 
Heyl |37] calculate linear modes of steady burning 
models, and find fuel is burned before the run- 
away. Bildsten has stressed that solving this puz- 
zle may require going beyond spherical symmetry, 
understanding both the burning front propaga- 
tion 21 J and the fuel covering fraction ^ . 

Some relevant observational facts are (i) regu- 
lar bursting seen in GS 1826 is well-understood as 
mixed H/He bursts, arguing that other bursters 
such as KS 1731-260 are in this burning regime 
when regular bursting is seen — if so (unless 
changes in covering fraction reverse the local ac- 



*Note that classifying short bursts as "helium" and long 
bursts as "hydrogen" is appropriate for bursts with simi- 
lar ignition column depths, but care should be taken when 
looking at bursts from a wide range of accretion rates. 
For example, at low accretion rates, the opposite behav- 
ior may be true, with long, energetic pure helium bursts 
and short hydrogen-triggered bursts. The best quantity 
to distinguish the composition is a. The lightcurves of a 
long mixed H/He burst and long pure He burst will also 
look very different; e.g. the latter will most likely show a 
large photospheric radius expansion episode. 
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cretion rate trend [3|), as M drops to the point 
where the recurrence time is or longer, a tran- 
sition at low M into pure helium bursts should 
be observed; (ii) large values a ~ 1000 are of- 
ten seen at high M, implying that both H and 
He are burning outside bursts; (iii) 4U 1820-30 
has a similar luminosity at which regular burst- 
ing stops — if the same mechanism applies, it 
must be composition- independent. 

Other sources show different behavior. In par- 
ticular, SAS-3 2 and later RXTE , .20.54) obser- 
vations show that bursts from 4U 1728-34 have 
helium like characteristics (a « 110 and dura- 
tions ~ 10 s), but become less energetic, stop 
showing PRE, and occur more often as accre- 
tion rate increases. The burst fluence is propor- 
tional to the recurrence time; very different to 
the « constant fluence expected for mixed H/He 
burning. Detailed models of the bursts from this 
source have not been made, but the observa- 
tions perhaps suggest accretion of He-rich mate- 
rial, with very little H burning between bursts. 

The difference in bursting behavior is corre- 
lated with the spin frequency measured during 
Type I bursts. KS 1731-260 has i/^pin = 
524 Hz; 4U 1728-34 has i^spin = 363 Hz. 

For both sources, burst oscillations are seen in 
the high accretion rate state only: for KS1731- 
260, this is mostly for PRE bursts; for 4U 1728- 
34, this is mostly for bursts without PRE. Muno 
and coworkers have shown that the correlation of 
the presence of burst oscillations with or without 
PRE applies to other sources with spins ~ 600 Hz 
or 300 Hz respectively [SEl- This suggests that 
a detailed analysis of burst properties would show 
that bursters with « 300 Hz spins are similar to 
4U 1728-34, and those with « 600 Hz are sim- 
ilar to KS 1731-260. It is not clear what this is 
telling us, perhaps that the two groups of systems 
have different evolutionary histories, and hence 
accreted compositions and burst properties. 

4. SUPERBURSTS 

4.1. Superbursts as carbon flashes 

A summary of observations of superbursts is 
given by Kuulkers, this volume (see also refs. |30| . 
|48| ) . Basic energetic arguments point to igni- 



tion of a layer of fuel much more massive than 
a typical Type I burst. For a nuclear energy re- 
lease of 1 MeV per nucleon, 10^^ g is needed 
to match the observed superburst energy of ~ 
10*^ ergs. At an accretion rate ~ 0.1 AfEdd 
appropriate for these sources, this mass is ac- 
creted in a few months. Since helium ignites at 
much lower masses at these M's (and indeed these 
sources show ordinary Type I bursts), the super- 
bursts require a new fuel. Theoretical attention 
has focused on carbon, which calculations show is 
present in the ashes of H/He burning, e.g. |8I46| . 

Carbon ignition on accreting neutron stars was 
studied by several authors in the 1970's 61 50 33^ 
(originally as a gamma-ray burst model), and 
more recently for rapidly accreting neutron stars 
(M > MEdd) by Brown and Bildsten 0. These 
calculations found ignition masses of ~ 10^^ g, 
giving characteristic energies lO^'^-lO^'* ergs, and 
recurrence times 10-100 years at 0.1 Medd- 
Therefore carbon flashes initially appeared too 
energetic to match the ^ 10^^ erg superbursts 
and probably too rare to match the number dis- 
covered (the superburst recurrence time is not 
well-constrained, but was at most « 6 years in 
one source 59 ). 

However, two pieces of physics act to re- 
duce the observed energies and recurrence times. 
Strohmayer & Brown 49 modeled the super- 
burst in 4U 1820-30 as ignition of a w 10^^ g 
layer with 30% carbon and 70% iron by mass. 
Because this source accretes helium-rich matter, 
stable burning is able to produce large amounts 
of carbon At the peak temperature of the 
flash (> 5 X 10^ K), neutrino emission is very 
efficient and Strohmayer & Brown found that 
~ 90% of the energy was lost as neutrinos, leav- 
ing « 10^^ erg to be radiated from the surface in 
the first few hours. 

The other superburst sources accrete hydrogen 
and helium, in which case protons readily cap- 
ture on carbon during hydrogen/helium burning. 
The amount of carbon produced is small, ^ 10% 
for stable burning, or < 1% for unstable burning 
1^ . However, Gumming & Bildsten ^21 showed 
that even a small amount of carbon is enough to 
power a superburst. The key physics is that the 
heavy elements that make up most of the mix- 



8 



hi 



10= 





1 1 ' 

\ 


1 ' 






\ 












\ 









10" 10» IQi" 10" 1012 1013 
Column Depth (g cm^^) 



Figure 5. Temperature profiles in the car- 
bon/heavy element layer (solid lines) for different 
carbon mass fractions. The short and long dashed 
lines show carbon ignition curves for Xc = 1-0 
and 0.1; the dotted line shows the pure helium 
triple alpha ignition curve for comparison. 



ture have a low thermal conductivity (this is es- 
pecially true for the very heavy nuclei made in 
the rp-process). This makes a steeper tempera- 
ture gradient, leading to ignition at a factor « 10 
lower in column depth. This is illustrated in Fig- 
ureEl which shows temperature profiles for differ- 
ent carbon mass fractions^. For Xc ~ 0.1, almost 
all of the w 10"*^ erg released leaves from the sur- 
face, since neutrino emission is not important at 
the lower peak temperature. The recurrence time 
is « 3 years for y w 10^^ g cm^^. 

Superburst sources have a small range of Af, 
K, 0.1-0.25 MEdd (Kuulkers, this volume). Gum- 
ming & Bildsten [HI showed that for M < 
0.1 A^Edd, the carbon burns stably before reach- 

take M = 0.1 MEddi a-^d assume the heating of the layer 
is by a flux from below of 0.1 MeV per nucleon. Brown 
171 has recently emphasised that the flux from below is 
very sensitive to the thermal conductivity and superfluid 
properties of the crust so that (unless the carbon layer 
has an internal heat source) this may be a promising way 
to constrain the thermal state of the interior. 
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Figure 6. Surface flux against time for cooling 
models for superbursts (Gumming & Macbeth, in 
preparation). Xc is the carbon mass fraction, y 
is the layer column depth (g cm~^). 



ing the ignition curve, explaining the lack of su- 
perbursts at low rates. Another possibility is that 
this is related directly to the change in Type I 
burst behavior observed at Lx ~ 2 x lO'^^ erg s~^ 
in bursters such as KS 1731-260 (§3.3; [HI). Be- 
low the transition, where the accreted fuel burns 
in regular Type I bursts, the carbon yield is likely 
very low; above, where the bursts are short and 
irregular and not burning all the fuel, stable burn- 
ing could produce significant amounts of carbon. 
An observational test is to compare bursts in 
sources with and without superbursts. The regu- 
lar burster GS 1826-24 should not show a super- 
burst in this picture. For M > 0.3 AfEdd, carbon 
flashes should still occur, but will be less ener- 
getic, and show less contrast with the accretion 
luminosity ^7]. Long duration « 10*^ erg bursts 
are seen from the high M source GX 17-1-2 j31j . 
but may be too frequent to be carbon flashes |E! • 

4.2. Time-evolution of the superburst 

A simple model of the superburst lightcurve 
is to assume the carbon burns instantaneously, 
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without significant vertical mixing. This may be 
a reasonable assumption since the burning occurs 
much more rapidly than the convective turnover 
time (opposite to normal Type I bursts). The 
thermal evolution of the layer can then be fol- 
lowed as it cools, by evolving the entropy equation 
in time (Gumming & Macbeth, in preparation). 
Figure El shows a series of lightcurves for differ- 
ent layer thickness (column depth y) and carbon 
fractions Xc^. 

At early times, as the outer parts of the 
layer thermally adjust, the radiative flux depends 
mostly on carbon fraction. At late times, after 
the cooling wave reaches the base of the layer, 
the flux falls off more steeply. A toy analytic 
model of the late time cooling is to take a slab 
with constant thermal diffusivity k, and perturb 
the temperature close to the surface, for example 
at a depth a. For a delta-function perturbation 
initially, the temperature evolution is given by 
the Green's function for this case, and the surface 
flux is F oc exp(— r/i)/t'^/^, where r is the ther- 
mal time at the initial heating depth r « o^/k. 
At late times, the flux decays as a power law 
oc close to the value found numerically. 

For y in the range 10^^-10^'^ g cm~^ and Xc — 
0.1-0.3, we find that ~ 10^^ ergs is radiated from 
the surface in the first several hours. The ob- 
served lightcurves agree well with these cooling 
models; a more detailed comparison is in progress 
to see whether Xq and y can be constrained. 

The calculations shown in Figure|Slassume that 
carbon burns to iron group after ignition, giving 
an energy release of w 1 MeV per nucleon. How- 
ever, Schatz et al. 01] showed that additional 
energy is released by photodisintegration of el- 
ements heavier than iron into iron group nuclei. 
The binding energy released in going from these 
heavy nuclei into iron is « 0.1 MeV per nucleon, 
ten times less than carbon burning into iron, but 
for Xc ^ 10%, this contribution dominates the 
energetics. This is an interesting additional way 
in which the nuclei made in the rp process directly 



°In reality, y and Xc are related by the ignition condi- 
tions, but for simplicity we treat y and Xc as free param- 
eters — this also gives a model independent approach in 
which Xc is just a measure of the deposited energy. For 
a similar approach see 1181 . 



affect the superburst energetics. 

4.3. Interaction betv^reen normal Type I 
bursts and superbursts 

Type I bursts disappear (are "quenched" ) for 
« weeks following the superburst. Gumming & 
Bildsten IJ, proposed that the cooling flux from 
the superburst stabilizes the hydrogen and helium 
burning. Assuming an exponential cooling law, 
they estimated a quenching time of ~ S^cooi ~ 
days, where icooi is the cooling time for the 
layer. This is a little shorter than the observed 
timescales of « weeks (Kuulkers, this volume). 
However, the cooling models in Figure El show a 
slow power law decay, F oc t~^/^, at late times 
rather than exponential. Inspection of Figure El 
shows that the flux reaches « lO^"* erg cm~^ s~^ 
in a time « 10 h (y/10^^ g cm~^), followed by 
a t^^/^ decay. Assuming the stabilizing flux is 
« 10^^ erg cm^^ s^^ |38l3j gives a quenching 
time 

iquonch ~ 22 tcooi ~ 9 days ( ^ ^ ) , (3) 

\ 10^^ g cm ^ / 

consistent with observed timescales. Turning this 
around, a measurement of the quenching time al- 
lows the thickness of the layer to be determined. 

Whenever the rise of the superburst has been 
seen, there is a precursor event that resembles a 
normal Type I burst. A likely explanation for 
this is that flux from deep carbon burning heats 
and ignites the layer of H/He accumulating on the 
surface. This scenario has yet to be investigated 
in detail, and in particular outstanding questions 
are the relative timing between precursor and su- 
perburst, and whether the vertical propagation 
of the burning during the rise triggers the precur- 
sor. However, simple predictions of this model 
are: (i) most superbursts should have a precursor 
(but not all, since a small fraction of the time the 
helium layer will be too thin to unstably ignite), 
(ii) a thicker helium layer will ignite with a lower 
flux, perhaps giving a correlation between pre- 
cursor energy and the time until the superburst 
rise, and (iii) the precursor should be weaker than 
Type I bursts at the same M, since the mass of 
the H/He layer is less than the critical mass. 

The final point is that simultaneous modelling 
of normal Type I bursts and superbursts gives 
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an important additional constraint on the mod- 
els of both J, 5 . For example, Strohmayer and 
Brown found a recurrence time of 13 years in their 
models of the superburst from 4U 1820-30 in 
which they inferred the accretion rate from the X- 
ray luminosity. However, using the regular Type 
I bursting to calibrate the accretion rate gives a 
shorter recurrence time, 1-2 years for pure he- 
lium, 5-10 years for 10% hydrogen [T3] . 

5. CONCLUSIONS 

The basic theory of nuclear burning on accret- 
ing neutron stars (first outlined over twenty years 
ago, e.g. has mixed success in explaining 

observed Type I X-ray burst properties. The 
extremely regular burster GS 1826-24 shows re- 
markable agreement with models of mixed H/He 
bursts. The observed variation of a with M 
strongly suggests that the accreted material has 
solar metallicity ,24 . Regular burst sequences 
from 4U 1820-30 agree well with accretion and 
burning of helium-rich matter. For this system, 
observations of burst fluence variations with re- 
currence time would constrain the amount of hy- 
drogen present in the accreted material |15| . rel- 
evant for testing evolutionary models. 

Superbursts fit naturally into the picture as ig- 
nition of fuel at a much larger depth than nor- 
mal Type I bursts. Calculations of H/He burn- 
ing show that the heavy ashes contain a small 
amount of carbon (~ 1-10%). Ignition calcula- 
tions show that this mixture will unstably ignite 
giving rise to a burst with properties matching 
observed properties of superbursts. Time depen- 
dent models of the superburst lightcurve show 
good agreement with the observed exponential 
decay times and energies. Because stable burn- 
ing produces more carbon than unstable burning 
|4f)| . a self-consistent model of carbon produc- 
tion and ignition may require understanding the 
nature of H/He burning at high accretion rates 
where bursts are irregular and fuel burns between 
bursts. There is much to be learned from the in- 
teractions between superbursts and normal Type 
I bursts. The properties of precursors and the 
time of cessation of bursting after a superburst 
will tell us much about the underlying layers. 



The global Type I bursting behavior remains 
the biggest puzzle, and may well need new physics 
to explain it. Another intriguing problem is 
how to explain bursts with ten minute recur- 
rence times (e.g. see EHI): not predicted by the- 
ory. There are many observational constraints 
which will allow progress to be made: variations 
of burst properties with accretion rate, shape 
of burst lightcurves, interactions between normal 
Type I bursts and superbursts, and the properties 
of burst oscillations. We are slowly learning more 
about bursts at higher or lower luminosities than 
traditional bursters (e.g. jl3l31| ). Understand- 
ing the rich phenomenology of bursts promises to 
teach us much about these neutron stars and the 
physical processes in their outer layers. 
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